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Nuclear Gluon Shadowing via Continuum Lepton Pairs in p + A at y/s = 200 AGeV 
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The dilepton spectra from open charm decay in p + A reactions at yfs = 200 AGeV are proposed to 
constrain gluon shadowing in nuclei. We show that the required measurements are feasible at RHIC. 



(N 
> 

O 
O 



X 



The spectrum of dileptons produced in heavy ion col- 
lisions at high energies has been proposed in the past 
to provide information about the dynamical evolution of 
quark-gluon plasmas (see, e.g., the review by Ruuskan- 
nen in ref. j^]). Unlike hadronic probes, that spectrum is 
sensitive to the earliest moments in the evolution, when 
the energy densities are an order of magnitude above 
the QCD confinement scale (1 GeV/fm'^). However, the 
small production cross sections for pairs with invariant 
mass above 1 GeV together with the large combinatorial 
background from decaying hadrons necessitate elaborate 
procedures to uncover the signal from the noise. The 
PHENIX detector , now under construction at the Rel- 
ativistic Heavy ion Collider (RHIC), is designed to mea- 
sure ee, e/z, and /z^ pairs to carry out this task. One of 
the important sources of background in the few GeV mass 
range arises from semileptonic decay of charmed hadrons 
{DD). As shown recently in ref. the expected ther- 
mal signals in that mass range may only reach 10% of the 
number of pairs from open charm decay. Therefore, spe- 
cial kinematical cuts and precise measurements will 
be needed to uncover the thermal signals 

There have been several attempts iQ,^ to take advan- 
tage of the large open charm background as a probe of the 
evolving gluon density. Most mid-rapidity charm pairs 
are produced via gluon fusion Therefore, the dilep- 
tons from charm decay carry information about the dis- 
tribution of primordial gluons before hadronization. In 
fact, the inside-outside cascade nature of such reactions 
greatly suppresses all sources of charm production 
except the initial perturbative QCD source. Therefore, 
the open charm background is dominated by the initial 
gluon fusion [gg cc) rate. 

The initial cc rate depends sensitively on the nuclear 
gluon structure function 0, gA{x,Q^). The quantity of 
fundamental interest [|l^,0 is the gluon shadowing func- 
tion 



Rg/Aix; Q^) = gA{x, Q'^)/AgN{x, Q^). 



(1) 



The point of this Letter, is to show that the A dependence 
of continuum dilepton pairs in the few GeV mass region 
provides a novel probe of that unknown gluon structure. 
We also show that the required measurements of ee,e/j-. 



and ^fx pair yields in p + A — > £^£~X at y/s = 200 
AGeV are not only experimentally feasible at RHIC but 
also that the open charm signal can be easily extracted 
via the proposed PHENIX detector. 

A key advantage of the continuum dilepton pairs from 
open charm decay over those from J/ip decay is that it is 
possible to test the applicability of the underlying QCD 
dynamics at a given fixed ^/s by checking for a particu- 
lar scaling property discussed below. In quarkonium pro- 
duction the mass is fixed and the required scaling can be 
checked only by varying the beam energy. In p+A —^ J/ip 
production the required scaling was unfortunately found 
to be violated in the energy range 20 < < 40 AGeV 
]l2[ thus precluding a determination of Ra ■ Possible ex- 
planations for the breakdown of QCD scaling for J/i/j 
production and its anomalous negative xp behavior [plp^] 
include interaction of next-to-leading order quarkonium 
Fock state with nuclear matter ||l^ , nuclear and comover 
Jjjp dissociation |lj] , and parton energy loss mechanisms 
|15| . It remains an open question whether the required 
scaling will set in at RHIC energies 60 < y/s < 200 
AGeV. 

Shadowing of the quark nuclear structure functions, 
qA{x,Q^), is well established from deep inelastic £A 
£X reactions [|l6|. For heavy nuclei, the quark structure 
functions are shadowed by a factor, Rq/^ix ^ 0.1, Q^) ~ 
1.1 — 0.1 A^/^, nearly independent of Q^. Perturbative 
QCD analysis |l^,|lll predicts that the gluon structure 
will also be shadowed at a; < 0.01 due to gluon recombi- 
nation processes. Therefore, a systematic measurement 
of gluon shadowing is of fundamental interest in under- 
standing the parton structure of nuclei. The gluon nu- 
clear structure is also of central importance in the field 
of nuclear collision since it controls the rate of mini-jet 
production that determines the total entropy produced 
at RHIC and higher energies |l^. We focus here on p+A 
collision rather than A-\-A\,o minimize the combinatorial 
TT, K decay backgrounds and other final state interaction 
effects. 

To demonstrate the sensitivity of open charm dilep- 
tons to gluon shadowing Rg/A{x,Q^), we calculate the 
pair spectrum from open charm decay in p + Au reactions 
at 200 GeV/A using for illustration two different gluon 
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shadowing functions. The first assumes that gluon shad- 
owing is identical to the measured quark shadowing, i.e. 
Rg/A{x, Q"^) = Rq/A{x, Q"^ = A GeV^), and is thus essen- 
tiaUy independent. This is the defauh assumption in- 
corporated in the HIJING Monte Carlo model . The 
second is taken from Eskola where i?g/^(a;,Q^) is 
computed using the modified GLAP evolution |10| start- 
ing from an assumed Rg/A consistent with the measured 
quark shadowing function at Q"^ = 4 GeV^. This second 
shadowing function depends on scale and differs for q and 
g. The two shadowing functions are shown in Fig. 1. 

The initial charm pair distribution from B + A nuclear 
collisions at impact parameter b is computed as in |^] : 



dp\dyidy2 



{Xa 7 Q 7 ) 



b,a 



di 



(2) 



where Ta/Aix,Q'^ ,r) = TA{r)fa/N{x,Q'^)Ra/A{x,Q'^,r) 
is the nuclear parton density function in terms 
of the known nucleon parton structure functions 
fa/N{x,Q'^), the nuclear thickness function ^^(r) — 
J dzuAiV + and the unknown impact parameter 
dependent shadowing function _Ra/^(a;, Q^, r). The con- 
ventional kinematic variables Xa,b sue the incoming par- 
ton light cone momentum fractions, t = —{pb —pi)'^, and 
the final c, c have rapidities yi , 2/2 and transverse mo- 
menta pi± = —p2± = P± ■ Next-to- leading order correc- 
tions to the lowest order parton cross sections daab / dt are 
approximately taken into account by a constant K factor 
We use the recent MRSA H parton structure 
functions. From a fit to low energy open charm produc- 
tion data in pp collisions, we fix rric = 1-4 GeV, if = 3, 
= s/2 for gg cc and s for qq cc. This choice of 
parameters leads to a charm pair cross section of 340/Li6 
for pp collisions at RHIC. For impact parameter averaged 
collisions, the integral over the transverse vector in eq.(|2) 
leads to a factor: BARb/Bixb,Q'^)Ra/A{xa,Q'^)/'jf^, 
where cr^"* is the inelastic B + A cross section and 
Ra/Aixa,Q^) is the impact parameter averaged shadow- 
ing function. 

In order to compute the DD pair distribution, a 
hadronization scheme for c D must be adopted. In 
e~^e~ cc DDX hadronization can be modelled via 
string fragmentation or fitted for example via the Peter- 
son fragmentation function (see (2^). The final D car- 
ries typically only a fraction ~ 0.7 of the original c mo- 
mentum. However, in pp collisions charm hadronization 
is complicated by the high density of partons produced 
during beam jet fragmentation. In this system recombi- 
nation or coalescence of the heavy quark with comoving 
partons provides another mechanism which in fact seems 
to be dominant at least at present energies. The inclusive 
pp DX data is best reproduced [plpj] with the delta 



function fragmentation D{z) = 5{1 — z) in all observed 
Xf and moderate pt regions. This observation can be 
understood in terms of a coalescence model if the coa- 
lescence radius is Pc ~ 400 MeV. We assume that hard 
fragmentation continues to be the dominant mechanism 
at RHIC energies, and hence no additional A dependent 
effects due to hadronization arise. However, this assump- 
tion must be tested experimentally. Our main dynamical 
assumption therefore is that the pp ^ D transverse mo- 
mentum distributions can be accurately reproduced from 
the QCD level rates assuming hard fragmentation as at 
present energies . At RHIC this can be checked either 
via single inclusive leptons or directly via Kir . 

With this assumption, the impact parameter averaged 
DD pair distribution in p + A is given by 



dNP^ 



KA E3E4 
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dp^dyzdy^ 

X! ^bfb/N{xb: Q'^)Xafa/N{Xa, Q'^)Ra/AiXa, -^^^^^ 

b,a 



dt 



(3) 



From this distribution of charmed mesons, we computed 
the dilepton spectrum via Monte Carlo using JETSET7.3 
|6| to decay the DD. 

While the lepton pair mass M and rapidity y fluctu- 
ate around a mean value for any fixed Mcs and j/ce, on 
the average (y) w (j/ce), and (Mcs) can be well approx- 
imated by a linear relation (Mcc) ~ PM + AIq over the 
lepton pair mass range 1 < M < 4GeV. For the delta 
function fragmentation case, f3 w 1.5 and Mq « 3.0 GeV 
are determined by the D-meson decay kinematics. 

Since gluon fusion dominates, there is an approximate 
scaling in terms of the light cone fraction xa that one of 
the gluon carries from nucleus A: 



lux A = -ycc + ln(Mcc/\/s) 

« -y + \n[il3AI + AIo)/V^ 



(4) 



The ratio of the dilepton dN/dMdy spectra in pA to 
scaled pp for different pair masses is thus expected to 
scale approximately as 



RP;^{M, y^~lnxA + ln[(/3M + Mo)/V^s] ) 
Rg/A{xA,Q^~il3M + Mof/2) 



1 dNP^ 
V dNll 



(5) 



where v = Aa^^/a^^ ^ A^/'^. In order to test gluon 
dominance and the accuracy of the above approximate 
scaling we compare in Fig. 1 the gluon shadowing func- 
tion to the above dilepton ratio from the Monte Carlo 
calculation. The solid curves are the input gluon shad- 
owings as a function of the Bjorken x. The other six 
curves are ratios of dilepton dN / dMdy spectra of shad- 
owed p + Au over those from unshadowed p + Au as given 
by eq. (j^) as a function of the scaling variable xa for three 
different dilepton masses. 
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FIG. 1. The dilepton dN/dMdy ratio curves (shadowed 
over unshadowed) at different masses obtained via Monte 
Carlo calculation are compared with the shadowing curves. 
The upper solid curve is Eskola gluon shadowing for 
q2 ^ iQ GeV^, and the lower solid curve is from HIJING 
[ p^ . The ratio curves are plotted as a function of the scal- 
ing variable logj^Q a;A from eq.(^). The scaled dilepton ratios 
reflect closely the input shadowing functions. 

Overall, the scaled ratio of lepton pair spectra approx- 
imates the shadowing function remarkably well. We con- 
clude that this ratio can therefore serve to map out gluon 
shadowing in nuclei. Note that in the case of Eskola shad- 
owing, even the dependence of the shadowing function 
is visible through the rise of the ratio curves with M in 
the small x region. 

Finally, it is important to estimate the dilepton back- 
ground to see if the proposed signal is experimentally 
feasible. Therefore we also calculated the signal and 
background dileptons for the PHENIX detector [|| taking 
into account the detector geometry and specific kinemat- 
ical cuts. The electron background is mainly due to tt*^ 
Dalitz and photon conversions. The background muon 
arises from random decays of charged pions and kaons. 
The electrons from 7r° Dalitz decay can be suppressed by 
small angle cut on dielectrons and the background 
muons are mainly suppressed by reducing the free decay 
volume. For the detector geometry, the electron arms 
cover pseudo-rapidity range —0.35 < rje < 0.35, and az- 
imuthal angle range ±(22.5°, 112.5°). The muon arm 
covers pseudo-rapidity range 1.15 < rj^ < 2.44 and al- 
most full azimuthal angle. For the kinematical cut, we 
take i?e > 1 GeV, > 2 GeV, and require that the 
relative azimuthal angle of the lepton pair ipi+e- > 90° 
in order to improve the signal-to-background ratio. 

For this background study, we use central p + Au col- 
lision with HIJING shadowing. We calculate the ee, e/i, 
and /z^ signal and backgrounds coming into the detec- 



tor, where the background electrons and muons are gen- 
erated from HIJING Monte Carlo calculation. We find 
(see Fig. 2) that the signal-to-background ratio for ee is 
very large, thus the dielectron signal from open charm 
decay is the easiest to extract. That ratio falls to about 
2.5 for e/i, and about 1/4 for /i/i. The pair rapidity distri- 
butions of the signal and backgrounds for central p -\- Au 
collision are shown in Fig. 2. Due to the detector geom- 
etry, ee, e^L, and /i// spectra cover pair rapidity regions 
centered at about 0, 1, and 2 (One can also study Au + p 
collision and measure pair rapidity regions centered at 
— 1 and —2). Like-sign subtraction should signifi- 
cantly reduce the noise, especially in the fj,fi channel. We 
conclude that the proposed measurement is feasible. 
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FIG. 2. Charm signal and backgrounds coming into the 
detector are plotted as a function of the lepton pair rapidity 
y. The solid curves are the signals from open Charm decay 
— labeled CC signal. The label C refers to the lepton coming 
from open Charm decay, label D refers to the electron coming 
from Dalitz and photon conversion, and label R refers to the 
muon coming from Random decay of pions and kaons. In (a) , 
the dielectron backgrounds are the dashed curve DC and the 
dot-dashed curve DD; and both are scaled up by a factor of 
100. In (b), the opposite-sign e/i backgrounds are the dotted 
curve CR, the dashed curve DR, and the dot-dashed curve 
DC. In (c), the dimuon backgrounds are the dashed curve 
RR and the dot-dashed curve RC. 

Further details of the calculations will be reported 
elsewhere For example, the Cronin effect and en- 

ergy loss are estimated to lead to distortions up to 10% 
at RHIC energies. From studies of the nuclear de- 
pendence of transverse momentum of J jt]} production, 
the typical increment, bp\_iyA), due to multiple colli- 
sions is limited to ^ 0.34 GeV^ even for the heavi- 
est nuclei |30|| . This momentum spread is shared be- 
tween the c and c, and distributed approximately as 
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g{5pA.) = e-''Pi/^V7^A2, A2 ~ 0.17 GeV2. For en- 
ergy loss, we assume that the charm quark loses an en- 
ergy 5E in the lab frame, and the charm quark p± is 
reduced by (1 — e), where e = 5E /[m_\_sinh.{y -\- yo)]- 
Combining these two effects, we may write for the charm 
quark that the final = (p™* -|- 5p±_)[l — e). There- 
fore the D-meson spectrum F{p±) = d?N/dm\ becomes 
F'ipi.) = J F{p^')g{6p^)dSp±. Taking F(p^) (x 6""^^, 
a ~ 1.3 GeV~^ from HIJING, and expanding the convo- 
lution to lowest order, the relative change of the D-meson 
spectrum is given by 

F'{p^_)/F{pi_) ~ 1 - ap^e + a^AV4 

This is also the relative change of the cc pair spectrum 
when Mcc = 2m± in case of yi = j/2, Pj_i = P±2- For 
m± - 3GeV, SE ^ lOGeV, coshyo - 100, the rel- 
ative change in the pair spectrum is estimated to be 
F'{p^)/Fip^)c^l -0.1(e-''-l). 

In summary, we calculated the Icpton pair spectra from 
open charm decay in two different shadowing scenarios. 
By scaling the ratios for different mass ranges according 
to eq.(^, we showed that the dilepton rapidity depen- 
dence of those ratios on xa reproduces well the underly- 
ing gluon shadowing function defined in eq.(^. Finally 
we showed that the measurements required to extract the 
gluon shadowing are experimentally feasible at RHIC. 
The above analysis depends on the validity of our basic 
assumption regarding the hard fragmentation of charm 
quarks, and this can be checked explicitly via the single 
inclusive measurements of D production. 

We conclude by emphasizing the importance of deter- 
mining gluon shadowing in p -f ^ to fix theoretically the 
initial conditions in A + A. In A + A the open charm 
decay is regarded as an annoying large background that 
must be subtracted to uncover the thermal signal. In 
p + A that background becomes the signal needed to de- 
termine the incident gluon flux in A + A. The continuum 
charm dileptons in p + A at RHIC are likely to provide a 
unique source of information on the low xa nuclear gluon 
structure at least until HERA is capable of accelerating 
heavy nuclei. 
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